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ABSTRACT 

Using detailed mock galaxy redshift surveys (MGRSs) we investigate the abundance 
and radial distribution of satellite galaxies. The mock surveys are constructed using 
large numerical simulations and the conditional luminosity function (CLF), and are 
compared against data from the Two Degree Field Galaxy Redshift Survey (2dFGRS). 
We use Monte Carlo Markov Chains to explore the full posterior distribution of the 
CLF parameter space, and show that the average relation between light and mass is 
tightly constrained and in excellent agreement with our previous models and with that 
of Vale & Ostriker. The radial number density distribution of satellite galaxies in the 
2dFGRS reveals a pronounced absence of satellites at small projected separations from 
their host galaxies. This is (at least partially) due to the overlap and merging of galaxy 
images in the APM catalogue. Because of the resulting close-pair incompleteness we 
are unfortunately unable to put meaningful constraints on the radial distribution of 
satellite galaxies; the data are consistent with a radial number density distribution 
that follows that of the dark matter particles, but we can not rule out alternatives 
with a constant number density core. Marginalizing over the full CLF parameter space, 
we show that in a ACDM concordance cosmology the observed abundances of host 
and satellite galaxies in the 2dFGRS indicate a power spectrum normalization of 
as — 0.7. The same cosmology but with erg = 0.9 is unable to simultanously match 
the abundances of host and satellite galaxies. This confirms our previous conclusions 
based on the pairwise peculiar velocity dispersions and the group multipliticy function. 

Key words: galaxies: formation — galaxies: halos — galaxies: fundamental param- 
eters — dark matter — cosmological parameters — methods: statistical 



1 INTRODUCTION 

In the hierarchical formation scenario, satellite galaxies are 
associated with dark matter subhaloes, which, some time 
in the past, were accreted by their current host (or parent) 
halo. As satellites orbit their host galaxies they are subjected 
to various forces that try to dissolve them: dynamical fric- 
tion, tides from the central object(s) and impulsive collisions 
with other satellites. A detailed understanding of the abun- 
dance and radial distribution of satellite systems, therefore, 
provides important constraints on the outcome of these var- 
ious physical processes which are an essential ingredient of 
galaxy formation. 

Traditionally, satellite galaxies have mainly been used 
as kinematic tracers of the dark matter potential well of the 
parent halo (e.g., Carlberg et al. 1996; Little & Tremaine 
1987; Evans et al. 2000; Zaritsky et al. 1993, 1997; McKay 
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et al. 2002; Prada et al. 2003; Brainerd & Specian 2003; 
van den Bosch et al. 2004b). Since satellite galaxies are dis- 
tributed over the entire halo, they are ideally suited to mea- 
sure the total virial mass, more so than for example rotation 
curves, which only probe the potential out to a fraction of 
the virial radius. The distribution of satellite galaxies, how- 
ever, has received considerably less attention. In addition 
to a few studies aimed at testing the claim by Holmberg 
(1969) that the azimuthal distribution of satellite galaxies is 
anisotropic with respect to the orientation of the host galaxy, 
the so-called Holmberg effect (Zaritsky et al. 1997; Zarit- 
sky & Gonzales 1999; Sales & Lambas 2004), relatively few 
studies have focussed on the radial distribution of satellite 
galaxies (but see Lake & Tremaine 1980; Vader & Sandage 
1991; Lorrimer et al. 1994; Willman et al. 2004). 

The interest in satellite galaxies has recently increased 
considerably, largely due to the dramatic increase in com- 
puting power that has made it possible to resolve dark mat- 
ter subhaloes in cosmological numerical simulations (Tor- 
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men 1997; Ghigna et al. 1998; Klypin et al. 1999; Moore et 
al. 1999; Stoehr et al. 2002; Kravtsov et al. 2003; De Lucia et 
al. 2004; Diemand, Moore & Stadel 2004; Gill et al. 2004a,b; 
Weller, Ostriker & Bode 2004; Reed et al. 2004). This has re- 
sulted in a number of detailed studies of the abundances and 
both the spatial and the velocity distribution of subhaloes. 
Since satellite galaxies are thought to be associated with 
these subhaloes, this has opened the possibility to compare 
the statistical properties of satellite galaxies directly with 
those of dark matter subhaloes, resulting in two apparent 
inconsistencies. 

First of all, the radial distribution of dark matter sub- 
haloes is found to be spatially anti-biased with respect to 
the dark matter. However, the observed distribution of clus- 
ter galaxies seems to accurately follow a Navarro, Frenk & 
White (1997, hereafter NFW) distribution (Beers & Tonry 
1986; Carlberg, Yee & EUingson 1997; van der Marel et 
al. 2000; Lin, Mohr & Stanford 2004). Although this might 
signal an inconsistency of the standard CDM framework, 
semi-analytical models of galaxy formation that take the 
evolution of dark matter subhaloes into account can repro- 
duce the observed distribution of cluster galaxies (Springel 
et al. 2001; Diaferio et al. 2001; Gao et al. 2004b). The rea- 
son is that galaxies are dense concentrations of baryons in 
the centers of dark matter (sub)-haloes, which are more re- 
silient to tidal disruption than (the outer parts of) their dark 
matter haloes. 

The second apparent inconsistency concerns the abun- 
dances of satellite galaxies. Although the subhalo mass func- 
tion, when normalized to the mass of the parent halo, is 
found to be virtually independent of halo mass (Moore et 
al. 1999; De Lucia et al. 2004; Diemand et al. 2004; Weller 
et al. 2004; but see Gao et al. 2004a), the satellite popu- 
lations of galaxy clusters are very different from those of 
galaxy sized haloes (Klypin et al. 1999; Moore et al. 1999; 
D'onghia & Lake 2003). This must be telling us something 
important about the physics of galaxy formation, and nu- 
merous studies have focussed on various mechanism to ex- 
plain this discrepancy between the number of predicted sub- 
haloes and observed satellites (e.g., Kauffmann, White & 
Guiderdoni 1993; BuUock, Kravtsov & Weinberg 2000; Ben- 
son et al. 2002) 

At this point in time there is a strong need for better 
observational constraints regarding the statistical properties 
of satellite galaxies. In particular, does the radial distribu- 
tion of satellite galaxies in galaxy sized haloes follow a NFW 
profile as for clusters, or does it reveal a constant density 
core, as found for dark matter subhaloes? How does the 
abundance of satellite galaxies depend on halo mass, or on 
the luminosity of the host galaxy? Is this mass-dependent 
abundance consistent with the full galaxy luminosity func- 
tion and the halo mass function? The conditional luminosity 
function (hereafter CLF) formalism developed by Yang, Mo 
& van den Bosch (2003) and van den Bosch, Yang &: Mo 
(2003a) is ideally suited to address these questions. It de- 
scribes how many galaxies of given luminosity reside in a 
halo of given mass. In this paper we therefore use the CLF, 
constrained using the abundances and clustering properties 
of galaxies in the Two-Degree Field Galaxy Redshift Sur- 
vey (2dFGRS, CoUess et al. 2001), to study the abundances 
and radial distribution of satellite galaxies. We construct de- 
tailed mock galaxy redshift surveys, based on the CLF, for 



direct comparison with the 2dFGRS data. We show that a 
close-pair incompleteness in the 2dFGRS prevents us from 
putting significant constraints on the radial number den- 
sity distribution of satellite galaxies, and that matching the 
abundances of host and satellite galaxies in the 2dFGRS re- 
quires a ACDM concordance cosmology with a relatively low 
value of the power-spectrum normalization parameter as- 

This paper is organized as follows. In Section 2 we de- 
scribe the CLF, and use a Monte-Carlo Markov Chain to 
fully sample the posterior distribution of the CLF. In Sec- 
tion 3 we describe how to use the CLF to construct detailed 
Mock Galaxy Redshift Surveys (MGRSs) for direct compar- 
ison with the 2dFGRS. The selection criteria for host and 
satellite galaxies are described in Section 4. Section 5 com- 
pares the projected, radial distribution of satellite galaxies 
extracted from the 2dFGRS with those obtained from our 
MGRSs, including a detailed discussion of incompleteness 
effects in the 2dFGRS. Section 6 compares the abundances 
of host and satellite galaxies between 2dFGRS and MGRS. 
We summarize our results in Section 7. 



2 THE CONDITIONAL LUMINOSITY 
FUNCTION 

Yang, Mo &: van den Bosch (2003) and van den Bosch, Yang 
& Mo (2003a) presented a new method to link the distribu- 
tion of galaxies to that of dark matter haloes. This method 
is based on modeling of the CLF, ${L\M)dL, which gives 
the average number of galaxies with luminosity L ± dL/2 
that reside in a halo of mass M. 

The CLF can be constrained by measurements of the 
abundances and clustering properties of galaxies. For in- 
stance, the CLF can be used to compute the luminosity 
function of galaxies 



$(L)dL : 



<&(L|M)n(M)dM 



(1) 



as well as the average bias of galaxies as function of lumi- 
nosity 



6(L) = 



1 



$(L) 



^{L\M)b{M)n{M)dM. 



(2) 



Here n{M) and fe(M) are the halo mass function and the 
halo bias, respectively. 

Throughout this paper we compute the halo mass func- 
tion using the form suggested by Sheth, Mo & Tormen 
(2001), which has been shown to be in excellent agreement 
with numerical simulations as long as halo masses are de- 
fined as the masses inside a sphere with an average overden- 
sity of about 180 (Jing 1998; Sheth & Tormen 1999; Jenkins 
et al. 2001; White 2002). Therefore, in what follows we con- 
sistently use that definition of halo mass when referring to 
M. The halo bias, b{M), is computed using the fitting for- 
mula of Seljak & Warren (2004). The linear power spectrum 
of density perturbations is computed using the transfer func- 
tion of Eisenstein & Hu (1998), which properly accounts for 
the baryons, while the evolved, non-linear power spectrum, 
required to compute the dark matter correlation function, is 
computed using the fitting formula of Smith et al. (2003). 

Throughout this paper we assume a flat ACDM 
cosmology with flm = 0.3, f^A = 0.7, h — 
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-ffo/(100 kms^^ Mpc^^) = 0.7 and with a scale-invariant 
initial power spectrum. Unless stated otherwise, we adopt a 
normalization of ag = 0.9. 



2.1 Parameterization 

The CLF is parameterized by a Schechter function: 

^ f L 

I* 



$(L|M)dL 



(f)' 



exp(— L/L*) AL, 



(3) 



where L* = L*(M), a = a{M) and !>* = l>*(M) are aU 
functions of halo mass M. We write the average, total mass- 
to-light ratio of a halo of mass M as 



{M/L), 



1 /M 

2 VT/o 



M 



+ 






(4) 



This parameterization has four free parameters: a character- 
istic mass Ml, for which the mass-to- light ratio is equal to 
(M/L)o, and two slopes, 71 and 72, that specify the behavior 
of {M/L) at the low and high mass ends, respectively. Moti- 
vated by observations (Bahcall, Lubin & Norman 1995; Bah- 
call et al. 2000; Sanderson & Ponman 2003), which indicate 
a flattening of {M/L)m on the scale of galaxy clusters, we 
set {M/L)m = (M/L)ci for haloes with M > 10"/i"^ M©. 

A similar parameterization is used for the characteristic 
luminosity L*{M): 



M 1 (M\ 

with 



M 



M \^3 



/(a) = 



r(a + 2) 
r(a + i,i) 



(5) 



(6) 



Here r(a;) is the Gamma function and r(a, x) the incomplete 
Gamma function. This parameterization has two additional 
free parameters: a characteristic mass M2 and a power-law 
slope 73. For a{M) we adopt a simple linear function of 
log(M), 



&{M) =Ql5+7) log(Ml5), 



(7) 



with M15 the halo mass in units of 10^^ /i^^M©, ais = 
<3(Afi5 = 1), and rj describes the change of the faint-end 
slope a with halo mass. Note that once a and L* are given, 
the normalization $* of the CLF is obtained through equa- 
tion (4), using the fact that the total (average) luminosity 
in a halo of mass M is given by 



{L)m= <l'(L|M)LdL = l'*Z*r(5 + 2). 



(8) 



Finally, we introduce the mass scale Mmin below which we 
set the CLF to zero; i.e., we assume that no stars form 
inside haloes with M < Mmin. Motivated by reionization 
considerations (see Yang et al. 2003 for details) we adopt 
Mmin = I0^h~^ Mq throughout. 



2.2 Parameter Fitting 

The CLF, as specified above, has a total of 8 free parame- 
ters: two characteristic masses; Mi and M2, three parame- 
ters that describe the various mass-dependencies 71 , 73 and 
77, two normalization for the mass-to-light ratio, (M/ L)o and 



{M/L)^i, and a normalization of the faint-end slope, 015. 
Note that 72 is not a free parameter as it derives from re- 
quiring continuity in {M/L)m across M = lO^''/i~"'^M0. The 
data that we use to constrain the CLF consists of the 2dF- 
GRS luminosity function of Madgwick et al. (2002) and the 
galaxy-galaxy correlation lengths as function of luminosity 
obtained from the 2dFGRS by Norberg et al. (2002a). The 
goodness-of-fit of each model is described by its x^ value 
defined by x^ = xl + Xn, "^i^h 



2 



and 



2 

Xra 






E 



$(LO-$(LO 



A$(LO 



ro{Li) ~ fo{Li) 



Afo(Li 



(9) 



(10) 



Here $ and fo are the observed quantities, and N,s> — 35 
and Nr = S are the number of data points for the LF and 
the correlation lengths, respectively. 

In Yang et al. (2003) and van den Bosch et al. (2003a) 
we presented a number of CLFs that accurately fit these 
data and that were based on different assumptions regard- 
ing the free parameters. Motivated by a number of indepen- 
dent observational constraints, the majority of these models 
were constrained to have {M/L)ci — 500h (M/L)0. How- 
ever, subsequent studies have shown that these CLF models 
predict too many rich galaxy groups (Yang et al. 2004b), 
and pairwise peculiar velocity dispersions that are too high 
(Yang et al. 2004a). Both these problems are alleviated 
by adopting a much higher cluster mass-to-light ratio of 
(M/L)ci ~ 900/1 (M/L) 0. 

In this paper we adopt a different approach. Rather 
than fixing (M/L)ci at a preferred value and using a min- 
imization routine to search our multi-dimensional parame- 
ter space for the best-fit parameters, we follow Yan, Madg- 
wick & White (2003) and use a Monte-Carlo Markov Chain 
(hereafter MCMC) to fully describe the likelihood function 
in our multi-dimensional parameter space. This will allow 
us to more accurately investigate the freedom in cluster 
mass-to-light ratio. Readers not familiar with, or interested 
in, Monte-Carlo Markov Chains are referred to Gamerman 
(1997) for details. 

We start our MCMC from model D in van den Bosch 
et al. (2004a) and allow a 'burn-in' of 10.000 random walk 
steps for the chain to equilibrate in the likelihood space. 
At any point in the chain we generate a new trial model 
by drawing the shifts in its eight free parameter from eight 
independent Gaussian distributions. The probability of ac- 
cepting the trial model is 



accept 



1.0 

exp[ 



<X 



2 
new 



ifXn 

Xoid)/2] ifx 



2 
ncvi' 



<Xold 
>Xold 



(11) 



with the x^ measures given by eq. (9) and (10). 

We construct a MCMC of 40 million steps, with an av- 
erage acceptance rate of ~ 12 percent. In order to suppress 
the correlation power between neighboring models in the 
chain, we thin the chain by a factor 20.000. This results in 
a final MCMC consisting of 2000 independent models that 
properly sample the full posterior distribution. The hatched 
histograms in Fig. 1 plot the resulting distributions of pa- 
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Figure 1. Constraints on the nine CLF parameters obtained from Monte-Carlo Markov Chains with 2000 independent samples. The 
hatched (non-hatched) histograms correspond to ACDM concordance cosmologies with erg = 0.9 (0.7), respectively. Vertical lines indicate 
the best-fit parameters for the erg = 0.9 cosmology. The median and 68 percent confidence intervals of these distributions are listed in 
Table 1. Masses and mass-to-light ratios are in units of h^^ TsAq and h (M/L)q, respectively. 



rameters, with the best-fit values indicated by a vertical line. 
The median and 68 percent confidence intervals of these dis- 
tributions are listed in Table 1 (model A0.9). Note that these 
parameters are similar to those of models A-D presented in 
van den Bosch et al. (2004a) '. 

In Fig. 2 we show scatter plots for some of the most 
tightly correlated pairs of parameters. The number density 
of points in these plots is directly proportional to the prob- 
ability density. The solid circles correspond to twelve spe- 
cific models to be discussed in more detail below. The most 
pronounced correlations are between {M/L)ci, {M/L)o and 
72: increasing the cluster mass-to- light ratio results in a de- 



T Small differences are due to the fact that here we use the halo 
bias fitting function of Seljak &; Warren (2004), rather than that 
of Sheth, Mo & Tormen (2001), and the transfer function of Eisen- 
stein & Hu (1998), rather than that of Efstathiou, Bond & White 
(1992). 



crease of the minimum mass-to- light ratio {M/L)o- This 
simply refiects the conservation of the number of galaxies 
(constrained by the LF): putting fewer galaxies in clusters 
requires an increase in the occupation number in lower mass 
haloes. Since 72 reflects the slope of {M/L)m in between Mi 
and the cluster scale, it is clear that 72 has to increase with 
increasing {M/L)ci and/or decreasing {M/L)o- Another im- 
portant correlation is between Mi and 71, which is a direct 
consequence of the shallow faint-end slope of the galaxy LF. 
Because of the much steeper low-mass end of the halo mass 
function, {M/L)m has to increase with decreasing mass. 
Lowering Mi results in a steeper increase, and therefore a 
larger value of 71 . 

Fig. 3 shows confidence levels on various quantities com- 
puted from our MCMC. The blue open circles with errorbars 
in the upper left and upper middle panels indicate the data 
used to constrain the models. The colored areas indicate the 
68 and 99 percent confidence levels on $(!/) and ro(-L) com- 
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Figure 2. Correlations between various CLF parameters of tlie 2000 samples in the erg = 0.9 MCMC (thin dots). The thick, solid dots 
correspond to the twelve CLF models listed in Table 1 which represent the best-fit models in twelve intervals of cluster mass-to-light 
ratio, (M/L)(.i. Masses and mass-to-light ratios are in units of h^^ Mq and h (M/L)q, respectively. 



puted from the MCMC. Note the good agreement with the 
data, indicating that the CLF can accurately match the ob- 
served abundances and clustering properties of galaxies in 
the 2dFGRS. We emphasize that this is not a trivial result, 
as the data can only be fitted for a certain combination of 
cosmological parameters (van den Bosch et al. 2003b). 

The upper right-hand panel of Fig. 3 plots the faint- 
end slope of the CLF as function of halo mass. The data 
clearly favors models in which a increases with decreasing 
halo mass. At around the cluster scale the models favor fairly 
steep faint-end slopes with —1.0 ^ 3 ^ — 1.5, which is in ex- 
cellent agreement with independent studies of the luminosity 
functions in a number of individual clusters (e.g., Sandage, 
Bingelli & Tammann 1985; Beijersbergen et al. 2002; Tren- 
tham & Hodgkin 2002; Trentham & TuUy 2002). The blue 
open circles with errorbars indicate a{M) for three different 
halo masses M obtained by Eke et al. (2004) from an anal- 



ysis of groups in the 2dFGRS, and are in good agreement 
with our model. 

The lower left-hand panel of Fig. 3 plots the relation 
between halo mass M and the total halo luminosity L, the 
expectation value of which follows from the CLF according 
to eq. (8). Note that the confidence levels are extremely 
tight, especially for the more massive haloes. The L{M) 
relation reveals a dramatic break at around M = M\ ~ 
7 X lQ^°h-^ Mq, with L oc M^'^^^-^ for M « Mi and 
L (X M°-^^*°°® for M » Ml (errorbars are obtained from 
the confidence levels on 71 and 72, respectively). In a recent 
study, Vale & Ostriker (2004), under the assumption of a 
monotonic L(M), obtained the relation between light and 
mass from a comparison of the galaxy luminosity function 
with the total halo mass function (counting both parent and 
subhaloes). They find that L{M) changes from L oc M^ to 
L oc M"'^ , in good agreement with our results. This is also 
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Figure 3. Posterior constraints on a number of quantities computed from the erg = 0.9 MCMC. The contours show the 68% and 99% 
confidence Umits from the marginaUzed distribution. Upper left-hand panel: The galaxy luminosity function; open circles with errorbars 
correspond to the 2dFGRS data from Madgwick et al. (2002). Upper middle panel: galaxy-galaxy correlation lengths as function of 
absolute magnitude; open circles with errorbars correspond to the 2dFGRS data from Norberg et al. (2002a). Upper right-hand panel: 
the faint-end slope of the CLF as function of halo mass; open circles with errorbars correspond to the 2dFGRS data from Eke et al. (2004). 
Lower left-hand panel: the total luminosity per halo as function of halo mass. The solid blue line corresponds to the model of Vale & 
Ostriker (2004), and is shown for comparison. Lower middle panel: the average mass-to-light ratio as function of halo mass. The various 
blue lines correspond to models A-D of van den Bosch et al. (2003a), as indicated. Lower right-hand panel: The average number of 
galaxies with Mjjj — 51og/i < —14 per halo as function of halo mass. See text for a detailed discussion. 



evident from a direct comparison of their results (blue solid 
line in the lower left panel of Fig. 3) with ours, which shows 
almost perfect agreement. It is extremely reassuring that 
two such wildly different methods yield results in such good 
agreement. This combined with the extremely tight confi- 
dence levels obtained from our CLF analysis suggests that 
we have established a remarkably robust connection between 
galaxy light and halo mass (at least for the concordance cos- 
mology adopted here). 

The lower, middle panel of Fig. 3 plots the correspond- 
ing mass-to-light ratios as function of halo mass. The pro- 
nounced minimum in AI/L indicates that galaxy forma- 
tion is most efficient in haloes with masses in the range 
5x 10"/i"^Mo < M < lO^^/i"^M0. For less massive haloes, 
the mass-to-light ratio increases drastically with decreasing 
halo mass, which is required in order to bring the steep slope 
of the halo mass function at low M in agreement with the 
relatively shallow faint-end slope of the observed LF. It in- 
dicates that galaxy formation needs to become extremely 
inefficient in haloes with M < 5 x 10^^'h~^ M© in order to 
prevent an overabundance of faint galaxies. The increase 



in {M/L)m from M ~ 10"/i"^ Mq to M ~ W^'^h'^ Mq 
is associated with the decreasing ability of the gas to cool 
with increasing halo mass (e.g., White & Rees 1978). Fi- 
nally, the sharp turn-over to a constant I^M/L)m for haloes 
with M > 10"/i"^ M© is a direct reflection of our CLF 
parameterization (see Section 2.1). The various lines corre- 
spond to models A-D from van den Bosch et al. (2003a), 
and are in excellent agreement with the confidence intervals 
obtained here. This indicates that the CLF models used in 
our previous work (Yang et al. 2004a,b; van den Bosch et 
al. 2003a,b; van den Bosch et al. 2004a,b; Mo et al. 2004; 
Wang et al. 2004) are perfectly consistent with the parame- 
ter constraints obtained here using the MCMC. 

Finally, the lower right-hand panel of Fig. 3 plots predic- 
tions for the number of galaxies (with Mbj — 51og/i < —14) as 
function of halo mass. These derive from the CLF according 



to 



{N)m = 



$(L|Af ) dL 



(12) 



with I/min the mininmm luminosity considered (Lmin = 
5.2 X W^h~^ L© in our case). Although the exact shape 
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and normalization of {N)m depends on Lmin (e.g., van den 
Bosch et al. 2003a), the power- law behavior combined with 
a shoulder plus break at low M is in good agreement with a 
number of studies based on halo occupation numbers (Sel- 
jak 2000; Scranton 2003; Berlind et al. 2003; Magliocchetti 
& Porciani 2003; Kravtsov et al. 2003). 



3 MOCK GALAXY REDSHIFT SURVEYS 

In order to properly interpret the 2dFGRS data on the 
abundance and radial distribution of satellite galaxies we 
construct detailed mock galaxy redshift surveys (hereafter 
MGRS). These have the advantage that (i) we know ex- 
actly the true abundances and the true radial distributions, 
(ii) we can model the effects of various observational biases, 
and (iii) we can use exactly the same host/satellite selection 
criteria as for the 2dFGRS, making the comparison with the 
data straightforward. 

To construct MGRSs two ingredients are required; a 
distribution of dark matter haloes and a description of how 
galaxies of different luminosity occupy haloes of different 
mass. For the former we use large numerical simulations, and 
for the latter the conditional luminosity function. Ideally we 
would construct a MGRS for each of the 2000 models in 
our MCMC. Since this is computationally too expensive we 
adopt an alternative method. We determine the minimum 
and maximum values of (M/L)ci in the entire MCMC, and 
split the interval in twelve equal size logarithmic bins. For 
each of these bins we then determine the model that yields 
the lowest x^ ■ The parameters of the resulting twelve mod- 
els are listed in Table 1, and are indicated by thick solid 
dots in Fig. 2. For each of these twelve CLFs we construct 
a MGRS, which we use to assess the uncertainties on the 
abundances and radial distribution of satellite galaxies due 
to the uncertainties in {M/L)ci- In what follows we refer to 
MGRS based on CLF n as Mn, where n is the ID listed in 
Column (1) of Table 1. 

Details about the construction of the MGRSs can 
be found in van den Bosch et al. (2004b) and Yang et 
al. (2004a). Here we briefly summarize the main ingredi- 
ents. The MGRSs are constructed by populating dark mat- 
ter haloes in a stack of large numerical simulations, kindly 
provided to us by Y.P. Jing (see Jing & Suto 2002). The 
number of galaxies per halo, and their luminosities, are ob- 
tained from the CLF. We assume that the brightest galaxy 
in each halo resides at rest at the halo center, while all other 
galaxies are in an isotropic steady-state equilibrium and fol- 
low a radial number density distribution given by 



"-'(^)°^(^)""(^+7fc)' 



(13) 



(limited to r < rvir) with r^ the characteristic radius of 
the NFW profile, and a and TZ two free parameters. Unless 
specifically stated otherwise we adopt a = TZ = 1 for which 
the number density distribution of satellite galaxies exactly 
follows the dark matter mass distribution. 

As described in detail in CoUess et al. (2001) and Nor- 
berg et al. (2002b), the parent catalogue of the 2dFGRS, 
which is extracted from the APM catalogue, is only 91% 
complete. In addition, because of problems with the tiling 
strategy and because of noisy spectra the level of complete- 



OS 




1.5 2 2.5 3 

log[0] (arcsec) 



3.5 



Figure 4. The solid line corresponds to the r(9) (see text for 
definition) obtained from the 2dFGRS by Hawkins et al. (2003). 
Note that r is significantly larger than unity for 9 < 100", indi- 
cating that the 2dFGRS has missed a significant fraction of close 
pairs (largely due to fiber collisions). The open circles with er- 
rorbars correspond to the r{6) obtained for MGRS M6, after we 
have removed close pairs. The good agreement with the results 
from the 2dFGRS indicates that we have constructed a MGRS 
with the same deficiency of close pairs as in the data. Note that 
this does not correct for a possible close-pair incompleteness in 
the parent catalogue of the 2dFGRS (see Section 5.2). 



ness of the 2dFGRS depends on position, apparent mag- 
nitude, and pair-separation. In the MGRSs used in Yang 
et al. (2004a), Wang et al. (2004) and van den Bosch et 
al. (2004b) we corrected for the average incompleteness of 
the 2dFGRS parent catalogue and for the position and mag- 
nitude dependent incompleteness, using detailed complete- 
ness maps provided by the 2dFGRS team, but no correction 
was made for the pair-separation dependency. For the pur- 
pose of investigating the radial distribution of satellite galax- 
ies, however, this close pair incompleteness needs to be prop- 
erly accounted for. Hawkins et al. (2003) computed the ratio 
r{8) = [1 + Wp{6)]/[1 + Wz{6)] between the angular correla- 
tion functions of the 2dFGRS parent catalogue, Wp{9), and 
that of the final redshift survey, w^ (9) . They found r to be 
significantly larger than unity, corresponding to a deficiency 
of pairs, for 9 < 100". In an attempt to model this close- 
pair deficiency in our MGRSs we proceed as follows. After 
having corrected for the position- and magnitude-dependent 
incompleteness, we compute the angular separations 9 be- 
tween all galaxy pairs and remove galaxies based on a prob- 
ability p{9), which we tune (by trial and error) so that we 
reproduce the r{9) obtained by Hawkins et al. (2003). As a 
last step we then remove a number of galaxies completely 
at random such that, together with the galaxies removed 
because of their angular separation to neighbors, we have 
removed 9 percent of all mock galaxies. This mimics the in- 
completeness level of the parent catalogue^ . The solid line in 



T Note that the sequential order in which the various complete- 
ness corrections have been applied is not entirely correct. How- 
ever, numerous tests have shown that this has an absolutely neg- 
ligible effect on our results. 
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Table 1. Conditional luminosity function parameters. 



ID 

(1) 



(M/LU 
(2) 



(M/L)o 
(3) 



logMi 
(4) 



logM2 
(5) 



71 
(6) 



72 
(7) 



73 
(8) 



(9) 



V 
(10) 



(11) 



1 


311.2 


167.1 


10.76 


12.45 


3.803 


0.176 


0.609 


-1.71 


-0.374 


70.53 


2 


353.2 


155.2 


10.61 


12.46 


4.310 


0.194 


0.650 


-1.62 


-0.314 


68.46 


3 


393.6 


135.8 


10.77 


12.22 


2.770 


0.236 


0.624 


-1.60 


-0.373 


65.97 


4 


441.6 


125.0 


10.72 


12.24 


3.141 


0.259 


0.646 


-1.55 


-0.330 


65.09 


5 


490.9 


115.1 


10.76 


12.10 


2.639 


0.287 


0.666 


-1.41 


-0.312 


64.34 


6 


508.2 


111.9 


10.77 


12.08 


2.594 


0.297 


0.654 


-1.46 


-0.333 


64.42 


7 


563.6 


102.1 


10.88 


11.93 


2.106 


0.334 


0.664 


-1.31 


-0.327 


65.18 


8 


632.4 


95.1 


10.93 


11.86 


1.970 


0.366 


0.665 


-1.29 


-0.342 


66.20 


9 


726.1 


85.7 


11.15 


11.66 


1.451 


0.431 


0.653 


-1.23 


-0.400 


67.95 


10 


798.0 


82.2 


11.21 


11.58 


1.331 


0.461 


0.639 


-1.25 


-0.445 


68.91 


11 


899.5 


76.4 


11.30 


11.64 


1.444 


0.507 


0.684 


-1.10 


-0.310 


72.01 


12 


1106.6 


71.4 


11.71 


11.48 


0.976 


0.652 


0.656 


-1.11 


-0.401 


72.99 



Ao.9 
Ao.7 



rnn + lSO 
O'J'J-100 



130 



1-15+23 



10 76+0-25 
iu.(0_g 21 



12 15+0-24 



2.92 



+2.11 
1.15 



29+°-''** 
^•^"-0.06 



0.66^ 



1 44+0.21 



-0.31 



+0.10 
0.10 



^""-100 



125 



+30 



10 71+0-21 



12.32^ 



3.46 



+2.18 
1.35 



0.24 



+0.08 
0.07 



0.62^ 



-1.64 



+0.21 



-0.35 



+0.10 



570+« 



98^ 



10.86«:«1 



11.96 



+0.10 



2.45 



+0.55 



0.34 



+0.02 
0.01 



0.63^ 



'1.35^ 



-0.35 



+0.08 



Parameters of CLF models. Column (1) lists the ID by which we refer to each CLF in the text. Columns (2) to (10) list the model 
parameters, and column (11) the value of x^ = X% + Xrn • Masses and mass-to-light ratios are in h~^M.Q and h (M/L)q , respectively. The 
first twelve lines (IDs 1-12) correspond to the best- fit models extracted from the erg = 0.9 MCMC for different bins in (M/L)^i. They 
are shown as thick solid dots in Fig. 2. The last three lines list the median and 68 percent confidence levels of the parameter probability 
distributions obtained from the MCMCs. Models Ao.g and Ao.7 correspond to the MCMC of the ACDM cosmologies with as = 0.9 and 
(Tg = 0.7, respectively. Model Ag*!^ is the same as Model A0.7, except that we have weighted the MCMC samples by exp(— ^2 72) (see 
Section 6.1). 



Fig. 4 shows r{9) for the 2dFGRS as obtained by Hawkins 
et al. (2003), while open circles with errorbars correspond to 
the r{9) obtained from a comparison of the angular correla- 
tion functions in MGRS M6 before and after the removal of 
close pairs. The good agreement with the 2dFGRS results 
demonstrates that we have managed to construct MGRSs 
with the same deficiency of close pairs as in the real 2dF- 
GRS. All MGRSs discussed in this paper have been cor- 
rected for this close-pair incompleteness. 



4 SELECTING HOST AND SATELLITE 
GALAXIES 

A galaxy is considered a potential host galaxy if it is at least 
fh times brighter than any other galaxy within a volume 
specified by Rp < Rh and \AV\ < {AV)h- Here Rp is the 
separation projected on the sky at the distance of the can- 
didate host, and AV is the line-of-sight velocity difference. 
Around each potential host galaxy, satellite galaxies are de- 
fined as those galaxies that are at least fs times fainter than 
their host and located within a volume with Rp < Rg and 
\AV\ < {AV)s- Host galaxies with zero satellite galaxies are 
removed from the list of hosts. 

In total, the selection of hosts and satellites thus de- 
pends on six free parameters: Rh, {AV)h and fh to spec- 
ify the population of host galaxies, and Rs, {AV)a and 
fs to specify the satellite galaxies. These parameters also 
determine the number of interlopers (defined as a galaxy 
not physically associated with the halo of the host galaxy) 
and non-central hosts (defined as a host galaxy that is not 
the brightest, central galaxy in its own halo). Minimizing 
the number of interlopers requires sufficiently small Rs and 



{AV)s- Minimizing the number of non-central hosts requires 
one to choose Rh, (AV')h and fh sufficiently large. Of course, 
each of these restrictions dramatically reduces the number 
of both hosts and satellites, making the statistical estimates 
more and more noisy. 

In van den Bosch et al. (2004b) we used our MGRS to 
optimize these selection criteria, aiming for large numbers of 
hosts and satellites (to increase signal-to- noise), a small frac- 
tion of interlopers, and a small fraction of non-central hosts. 
It was shown that an adaptive selection criterion, for which 
Rh, {AV)h and Rs are made dependent of the luminosity of 
the host-candidate, was most successful. Motivated by these 
findings we adopt: fh ~ fs = 1, {AV)h ~ 1000(J2oo kms^^, 
(AV)s = 2000 kms"\ Rh = 0.80-200/1"^ Mpc, and Rs = 
0.15a2oo/i"^ Mpc. Here a2oo is the satellite velocity disper- 
sion, (Tsat(ihost), in units of 200 kms^^, for which wo adopt 



logo-sat = logaio + ailogl/io + a2(logLi 



(14) 



Here Lio = Lhost/10 
-1 



h~ L© and 0-10 = o-sat(-Z^io)- We set 
aio = 200kms^^, ai = 0.5 and 02 = 0.1, which is close to the 
parameters obtained for the 2dFGRS by van den Bosch et 
al. (2004b). Note that details regarding the selection criteria 
of host and satellite galaxies are not very important as long 
as the data and model (MGRS) are treated in the same way. 
For the data we use the final, public data release from 
the 2dFGRS, restricting ourselves only to galaxies with red- 
shifts 0.01 < 2 < 0.15 in the North Galactic Pole and South 
Galactic Pole subsamples with a redshift quality parameter 
g > 3. This leaves a grand total of 146735 galaxies with a 
typical rms redshift error of 85 kms"^ (CoUess et al. 2001). 
Absolute magnitudes for galaxies in the 2dFGRS are com- 
puted using the K-corrections of Madgwick et al. (2002) . Ap- 
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Figure 5. The projected number density distributions of satellite galaxies as function of projected radius. Results are shown for three 
different bins in host-galaxy luminosity, indicated at the top of each column. In each panel, the hatched histogram corresponds to the 
2dFGRS data, while the various lines have been obtained from the MGRSs. The errorbars (only shown in the upper panels for clarity) 
indicate Poisson errors in the numbers of host and satellite galaxies. The upper panels compare the 2dFGRS data to three MGRSs (Ml, 
M6, and M12) that span the entire range in cluster mass-to-light ratio (see Table 1 for parameters). The middle row of panels compare 
the 2dFGRS data to MGRSs based on CLF 6, but in which the intrinsic, spatial distribution of satellite galaxies is modeled with a 
constant number density core (eq. [13] with = 0) with 7^ = 1 (M6a), 2 (M6b), and 3 (M6c). For comparison the fiducial MGRS M6, 
for which a = 7?. = 1 is also shown (solid lines). Finally, the lower panels compare the 2dFGRS data to MGRSs based on CLF 6, but 
in which we have corrected for the size-projection incompleteness using Rio = 0.010/i~ Mpc (M6d), Rio = 0.015/i~ Mpc (M6e), and 
RiO = 0.020h-i Mpc (M6f). See text for detailed discussion. 



plying our host-satellite selection criterion yields 8737 host 
galaxies and 13738 satellite galaxies. 



5 THE RADIAL DISTRIBUTION OF 
SATELLITE GALAXIES 

The hatched histograms in Fig. 5 show the projected, radial, 
number density distribution, Esat(-Rp) (in /i^/Mpc^/host) 
of 2dFGRS satellite galaxies around hosts in three different 
host magnitude bins. Note that the Ssat(fip) reveal both an 
inner and an outer break. The outer break radius increases 
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with increasing host luminosity, and is mainly due to our se- 
lection criteria. The inner break radius, however, indicates a 
pronounced absence of satellite galaxies at small projected 
radii from the host galaxy. In order to interpret these find- 
ings in terms of the actual three-dimensional number density 
distribution of satellite galaxies we compare these findings 
to those obtained from our MGRSs. 

The three non-hatched histograms in the upper panels 
of Fig. 5 depict the Ssat(-Rp) as obtained from MGRSs Ml 
(dotted lines), M6 (solid lines), and M12 (dashed lines). Note 
that all three MGRSs yield virtually identical radial num- 
ber density distributions. In fact, this holds for all twelve 
MGRSs listed in Table 1: for clarity, however, only three are 
shown. Although the projected number density distributions 
of the MGRSs nicely match the 2dFGRS data at large pro- 
jected radii, they severely overpredict T,sa.t{Rp) at small Rp. 
Rather than a pronounced inner break radius, Ssat(i?p) in- 
creases continuously down to the smallest radii shown. Espe- 
cially for the brightest host galaxies the discrepancy between 
MGRSs and 2dFGRS is dramatic. Either satellite galaxies 
are spatially anti-biased with respect to the dark matter, or 
the 2dFGRS has somehow missed a large number of satel- 
lite galaxies at small projected radii from their host galaxies. 
Below we investigate both options. 



5.1 Spatial anti-bias 

When constructing the MGRSs we have thus far assumed 
that satellite galaxies follow the same number density dis- 
tribution as dark matter particles (i.e., eq [13] with a = 
TZ = 1.0). The resulting overestimate of the projected num- 
ber density of satellite galaxies at small Rp, however, might 
indicate that a < 1.0 and/or that 7^ > 1. To test these ideas 
we construct three MGRSs using CLF 6. These MGRSs, 
termed M6a, M6b and M6c, all have a — 0.0, and only 
differ in the value of TZ, which is set to 1.0, 2.0, and 3.0, 
respectively. The projected number density distributions of 
their satellite galaxies are shown in the middle row of panels 
of Fig. 5. For comparison we also plot the T,Ba.t{Rp) of M6 
(solid lines). As expected, lowering a and increasing 7?. both 
result in a decrease of Ssat at small Rp. In fact, MGRS M6c 
fits most of the 2dFGRS data reasonably well, although it 
fails to reproduce the pronounced inner break evident in the 
data. 

Taking these results at face value seems to indicate that 
satellite galaxies are spatially anti-biased with respect to the 
dark matter mass distribution on small scales (20 to 30 per- 
cent of the virial radius). This is in wonderful agreement 
with high resolution numerical A'^-body simulations, which 
reveal a very similar spatial anti-bias for dark matter sub- 
haloes (Ghigna et al. 1998, 2000; Cohn et al. 1999; Okamoto 
& Habe 1999; Springel et al. 2001; De Lucia et al. 2004; Die- 
mand et al. 2004) , but inconsistent with observations of the 
number density distribution of galaxies in clusters (Beers & 
Tonry 1986; Carlberg et al. 1997; van der Marel et al. 2000; 
Diemand et al. 2004; Lin, Mohr & Stanford 2004) and with 
predictions of semi-analytical models of galaxy formation 
(Springel et al. 2001; Diaferio et al. 2001; Gao et al. 2004b). 



5.2 Observational projection effects 

This puzzling inconsistency calls for a closer examination 
of the selection effects in the 2dFGRS. Could the absence 
of satellite galaxies at small projected radii from their host 
galaxies be due to selection effects? Note that it can not be 
explained as due to the fiber-collisions. As shown in Fig. 4, 
and discussed in Section 3, these effects are accurately ac- 
counted for in our MGRSs. One effect that has not yet been 
accounted for, however, is related to the fact that galaxies 
have a finite size, which causes them to overlap in projec- 
tion. Whenever this happens only the brightest of the two 
will be recognized as a galaxy and included in the survey. 
Indeed, as shown by Cole et al. (2001) using a compari- 
son of the 2dFGRS with the Two Micron All Sky Survey 
(2MASS), about 4.5 percent of all galaxies is missed in the 
2dFGRS parent catalogue because of merged or close images 
(explaining about half of the incompleteness; the other half 
being due to incorrect star-galaxy separation) 

In order to take this size-projection effect into account 
we proceed as follows. We model the characteristic size of a 
galaxy as 



R. 



gal 



Ri 



L 



101«/i-2 Lq 



(15) 



and define the critical projection angle Smax ~ -Rgai/^A, 
with Da the angular distance of the galaxy. Based on the 
data on disk galaxies in Courteau et al. (2004) we adopt 
( — 1/3, but we note that this assumption is not very impor- 
tant. For example, we have also experimented with ( — 1/2 
finding very similar results as those described here. 

We construct three MGRSs in which we take the size- 
projection effect into account using Rio — O.OlO/i^^ Mpc 
(M6d), Rio = 0.015/1^^ Mpc (M6e), and Rw = 
0.020ft^^ Mpc (M6f). We follow the same procedure as de- 
scribed in Section 3, except that this time, after having cor- 
rected for the position- and magnitude-dependent incom- 
pleteness, we compute the angular separations 9 between all 
galaxy pairs and remove the faintest galaxy from all pairs 
for which 6 < ^max. Typically this removes about 2 to 4 
percent of all galaxies, which is comparable to the incom- 
pleteness level in the 2dFGRS parent catalogue due to this 
effect as found by Cole et al. (2001). Next we correct for 
the fiber collisions, using the same trial-and-error method 
as described in Section 3. Finally we remove a number of 
galaxies completely at random to bring the total fraction of 
removed galaxies to 9 percent. 

The lower panels Fig. 5 plot the Esat(i?p) of MGRSs 
M6d (dashed lines), M6e (dot-dashed lines), and M6f (dot- 
ted lines). Note that these MGRSs have a = TZ = 1.0, as 
for our fiducial MGRSs. For comparison, the solid lines plot 
the projected number density distributions of MGRS M6, 
for which no correction for the size-projection effect has 
been applied. Clearly, the size-projection effect as modeled 
here has a dramatic impact on Eaat(7?p) at small projected 
radii. In fact, it reproduces the inner cut-off observed for 
the 2dFGRS data, indicating that, as first noticed by Cole 
et al. (2001), the 2dFGRS lacks a significant fraction of close 
galaxy pairs due to projection effects. 

Unfortunately, unless we can model this particular in- 
completeness in more detail as done here (i.e., have indepen- 
dent constraints on Rio and ("), we can not use the 2dFGRS 
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Figure 6. Abundances of host and satellite galaxies. From left to right, the different panels in a given row plot the ratios A^ijost /^tot , 
Nsax/Ntot and Nsat / Ni^^st > a-H a-s function of cluster mass-to-light ratio {M/L)^i (in units of h (M/L)©). Gray, horizontal bars indicate 
the ratios obtained from the 2dFGRS. Upper panels: solid triangles and open circles indicate the ratios obtained from the MGRSs for 
the twelve CLFs listed in Table 1 without and with correction for the size-projection effects, respectively. Errorbars, computed assuming 
Poissonian errors on Ntot, Afjiost ^nd Afgat, are smaller than the symbols, and are therefore not shown. The solid line indicates the 
ratios obtained from the analytical estimates (eq. [16]-[21]) after taking account of the correction factors /host a^nd /sat (see text). The 
good agreement with the open circles indicates that these analytical estimates can be used to compute the abundances of host and 
satellite galaxies for any CLF model without the need to construct a detailed MGRS. Middle panels: the ratios Af^ost /^tot , Asat/A^tot 
and Afsat/A^host for all 2000 samples in the erg = 0.9 MCMC, obtained using the approximate, analytical method. Note that none of 
these models can simultaneously match all three ratios obtained from the 2dFGRS. Lower panels: Same as middle panels, but this time 
showing the results for the erg = 0.7 MCMC. Note that at around log[(M/L)ci] ~ 2.75 these model predictions simultaneously match all 
three 2dFGRS ratios. 



data to meaningfully constrain the actual radial density dis- 
tribution of satellite galaxies. Our models suggest that the 
data is consistent with no spatial (anti)-bias as long as the 
size-projection effects are taken into account, but we can not 
rule out that a — 0.0 and/or that TZ > 1.0. For example, 
models with both a = 0.0 and in which we correct for the 
size-projection effect are virtually indistinguisahble from the 



same models but with a = 1.0. It remains to be seen whether 
the higher spatial resolution of the SDSS will allow a more 
in-depth investigation of the radial distribution of satellite 
galaxies. 
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6 THE ABUNDANCE OF HOST AND 
SATELLITE GALAXIES 

We now focus on the abundances of host and satelhte galax- 
ies. For the 2dFGRS we obtain host and satellite fractions 
of AThost/A^tot = 0.060 (0.060, 0.059) and iV.at/iVhost = 0.094 
(0.096,0.092), where the numbers in brackets indicate the 
fractions obtained using only the NGP and SGP data, re- 
spectively. We compare these numbers with those obtained 
from the MGRSs. The solid triangles in the upper panels of 
Fig. 6 plot A''host/A''tot, A''sat/A^tot, and A'"sat/A''host as func- 
tion of (M/L)ci for the twelve MGRSs listed in Table 1. The 
grey, horizontal bars indicate the 2dFGRS results. Clearly, 
the fraction of satellite galaxies in the MGRS is much too 
high compared to the 2dFGRS, independent of {M/L)c\- 
The fraction of host galaxies is in reasonable agreement 
with the 2dFGRS (though somewhat too high) , but only for 
(M/L)ci < 650/i (M/L)0. Matching the number of satellites 
per host, however, favors much higher values for {M/L)ci- In 
short, none of the twelve MGRS can simultaneously match 
the abundances of host and satellite galaxies in the 2dFGRS. 

The open circles in the upper panels of Fig. 6 corre- 
spond to the same twelve MGRSs, but this time corrected 
for the size-projection effects discussed in Section 5.2 using 
Rw = 0.015/i"^ Mpc and C, = 1/3. This lowers the overall 
abundances of both hosts and (mainly) satellites, but still, 
none of the twelve MGRS can simultaneously match the 
abundances of host and satellite galaxies in the 2dFGRS. 
Although the host fraction is in good agreement with the 
data for (M/L)ci ^ 650/i (M/L)0, matching the number of 
satellites per host requires (M/L)ci > 800/i (M/L)©. 

In order to investigate this failure of the models in more 
detail, and to explore all freedom in the CLF parameters, 
ideally one would construct a MGRS for each of the 2000 
models in our MCMC. This, however, is not feasible com- 
putationally, and we therefore use an approximate, analyti- 
cal method instead. The total number of galaxies in a flux- 
limited survey is given by 



(ATs, 



A^t, 



AQ. 



dz / dMn{M){N)M,: 



(16) 



Here Q. is the solid angle of sky of the survey, AV is the 
differential volume element, Zmin and Zmax are the survey 
redshift limits (zmin = 0.01 and Zmax = 0.15 in our case), 
and (N)m,z is the average number of galaxies per halo of 
mass M at redshift z which follows from the CLF according 
to 



{N), 



A^) 



$(L|M)dL 



(17) 



Jz) 



with Lmin{z) and Lmax(z) the minimum and maximum lu- 
minosities of a galaxy at redshift z that makes the apparent 
magnitude limits of the survey (we adopt 15.0 < mtj < 19.3 
for the 2dFGRS). Similarly, the number of satellite galaxies 
follows from 



A^sat = 



with 



AQ 



dz / dMn(M){Ar,at)M,; 



(18) 



{N)m,z - 1 if {N)m,. > 1 
otherwise 



(19) 



The number of host galaxies in a flux-limited survey, finally, 
is given by 

A^host = AQ Az AMn{M)w{M,z) (20) 

where the weight function w{M, z) is given by 



w{M,z) = 



1 - CXp(-(A'sat>M,.) if {N)m,. > 1 







otherwise 



(21) 



This derives from the fact that only hosts with at least one 
satellite galaxy are counted as host galaxies, and using the 
fact that the number of satellite galaxies follows a Poissonian 
distribution (see van den Bosch et al. 2004b). 

Using these equations we compute the abundances of 
host and satellite galaxies for each of the 2000 CLF mod- 
els in the MCMC. However, this does not take care of in- 
completeness effects, fiber collisions, or the size-projection 
effect. In addition, this assumes that all satellites are se- 
lected with zero interlopers. Instead, the host-satellite se- 
lection criterion used for the 2dFGRS and the MGRSs is 
not 100 percent complete and yields about 15 percent in- 
terlopers (van den Bosch et al. 2004b). To correct for all 
these effects/shortcomings we proceed as follows. We mul- 
tiply the A'^host and A^sat computed using eq. (20) and (18) 
with correction factors /host and /sat, respectively, which 
we calibrate using the results of the 12 MGRSs (that have 
been corrected for the size-projection effects) shown in Fig. 6 
(open circles). We flnd a good match for /host ~ 0.8 and 
/sat = 0.57 + 0.25[log(M/L)ci - 3.0], the results of which 
are indicated by the solid line in the upper panels of Fig. 6. 
In what follows we assume that the same correction factors 
apply to all models in our MCMC. 

Although only approximate, eq. (16)-(21) combined 
with this simple scaling allows us to make predictions for 
the fractions of host and satellite galaxies for each of the 
2000 CLF models in our MGRSs. The results are shown 
in the middle row of panels of Fig. 6. The overall behav- 
ior is similar to that of the 12 MGRSs shown in the up- 
per panels: matching the fraction of host galaxies requires 
(M/L) ci <, 650h (M/L)0, whereas matching the number of 
satellites per host basically requires the opposite. Mean- 
while, the number of satellite galaxies is systematically too 
high. Although there are a few models that match the frac- 
tion of satellite galaxies of the 2dFGRS, as we show below, 
these models can not simultaneously match the fraction of 
host galaxies. Therefore, we conclude that the abundances 
of host and satellite galaxies in the 2dFGRS can not be 
reproduced within the ACDM concordance cosmology with 
0-8 = 0.9. 

6.1 Evidence for a low value of ag 

Using MGRSs similar to those presented here, Yang et 
al. (2004a) have shown that a ACDM concordance cosmol- 
ogy with (M/L)ci = 500h (M/L)© predicts too much clus- 
tering power on small scales, and pairwise peculiar velocity 
dispersions that are too high. In Yang et al. (2004b) it was 
shown that the same model predicts too many large groups 
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Figure 7. Same as Fig. 2 except tliat tliis time we plot parameter correlations for the MCMC constructed assuming a ACDM cosmology 
with (Tg = 0.7. The thick solid dots indicate those models for which x^b ^ ^' indicating a good match to the observed abundances of 
host and satellite galaxies. Note that these extra constraints severely restrict a number of the CLF parameters; cf., the parameters of 
models Aq.t and A^^^ in Table 1. 



of galaxies. All these problems indicate that there are too 
many galaxies in (massive) clusters compared to observa- 
tions. As demonstrated in Yang et al. (2004a, b), this can 
be remedied by adopting a high {M/L)ci, which results in 
fewer galaxies per cluster. This agrees perfectly with the 
fact that matching the number of satellites per host re- 
quires a similarly high (M/L)ci. However, the failure to si- 
multaneously match the separate abundances of host and 
satellite galaxies, and the independent observational con- 
straints on cluster mass-to-light ratios which indicate that 
(M/L)ci ~ (400 ± 100)h (M/L)0 (Carlberg et al. 1996; 
Fukugita, Hogan & Peebles 1998; BahcaU et al. 2000), signal 
a clear shortcoming of this model. 

An alternative solution to the apparent overabundance 
of cluster galaxies is to lower the abundance of cluster- 



sized haloes. This is most easily established by lowering 
the power-spectrum normalization ug. As shown by Yang 
et al. (2004a,b), a ACDM cosmology with erg = 0.7 and 
(A//L)ci = 500h (M/L)0 can equally well match the clus- 
tering data, the pairwise peculiar velocities, and the abun- 
dances of galaxy groups as a concordance cosmology with 
erg = 0.9 and {M/L)ci = 900h(M/L)Q. In addition, as shown 
in van den Bosch et al. (2003b), this cosmology is in good 
agreement with the WMAP data and may even alleviate 
some problems with the ACDM cosmology regarding the 
concentration of dark matter haloes. It is interesting, there- 
fore, to investigate whether such a model can also simulta- 
neously match the abundances of host and satellite galaxies 
in the 2dFGRS. 

In order to test this we construct a MCMC for the 
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ACDM cosmology but with as — 0.7. We use the same 
number of steps and the same thinning factor as for the 
MCMC described in Section 2.2. The resulting distribution 
of parameters is indicated by the non-hatched histograms 
in Fig. 1. The corresponding median and 68 percent con- 
fidence intervals are listed in Table 1 (model A0.7), Cor- 
relations amongst various parameters are shown in Fig. 7. 
Compared to the distributions for the concordance cosmol- 
ogy with (78 = 0.9 the main differences are a reduction of 
the mean {M/L)ci, 015 and 73. Other parameters, notably 
Ml and 71, are extremely insensitive to as- 

Unfortunately we do not have numerical simulations for 
a ACDM cosmology with o-g = 0.7, so that we can not con- 
struct MGRSs for this cosmology (but see Yang et al. 2004a). 
Nevertheless, we can use the scaling parameters /host and 
/sat, under the assumption that they are also valid for this 
cosmology, to predict the abundances of host and satellite 
galaxies for the low-erg cosmology using eq. (16)-(21). The 
resulting host-satellite fractions are shown in the lower pan- 
els of Fig. 6. Compared to the same results for the erg = 0.9 
cosmology (panels in middle row), the satellite fractions have 
been reduced, bringing them in better agreement with the 
2dFGRS results. In order to make the comparison between 
the two different ag models more quantitative, and to in- 
vestigate whether any model can simultaneously match the 
abundances of host and satellite galaxies, we introduce the 
goodness-of-fit measure 



2 _ 

Aab 



iVto 



0.060 



0.002 



+ 



Aftot 



0.094 



0.003 



(22) 



where the numbers in the numerators are the 2dFGRS ra- 
tios, and those in the denominators are the standard devi- 
ations due to cosmic variance that we obtain from a set of 
independent MGRSs using different simulation boxes (see 
Yang et al. 2004a) J. For the ag = 0.7 cosmology, the num- 
ber of samples (out of a total of 2000) with xlb < (2, 4, 6) is 
(11,34,57). For the as — 0.9 cosmology, these numbers are 
(0, 0, 0). Clearly, the erg = 0.7 models are far more succesful 
in sim,ultaneously fitting the abundances of host and satellite 
galaxies than those with erg = 0.9. 

The thick, solid dots in Fig. 7 indicate the 34 CLF mod- 
els in the erg = 0.7 cosmology for which Xab < 4. Note how 
they are clustered together in parameter spaces. This means 
that using the observed abundances of host and satellite 
galaxies in the 2dFGRS as extra constraints allows a sig- 
nificant tightening of the constraints on the CLF param- 
eters. We can take account of these additional constraints 
by weighting each sample in the MCMC with exp(— Xab/2). 
This yields the median and 68 percent confidence intervals 
listed under model Kq^-j in Table 1. Note that the constraints 
on some of the CLF parameters, notably {M/L)o, (M/L)ci 
and Ml, are now much tighter (cf. model A0.7). 

Although the above results strongly favor a ACDM cos- 
mology with relatively low power spectrum normalization, 
we caution that they are obtained assuming that the scal- 
ing parameters /host and /sat, determined for a small sub- 



3 Note that these are very similar to the standard deviations 
obtainded from a comparison of the host/sateUite fractions in the 
2dFGRS NGP and SGP 



set of all MCMC samples with erg = 0.9, are valid for all 
MCMC samples, even those with erg = 0.7. Testing the ac- 
curacy of this assumption requires similar A'^-body simula- 
tions as those used here, but for the ACDM cosmology with 
erg = 0.7, and the construction of a large number of differ- 
ent MGRSs to probe the full CLF parameter space. Given 
our limited computational resources, we are unfortunately 
unable to perform these tests. 



7 SUMMARY 

Comparing the radial distribution and abundances of dark 
matter subhaloes in simulations with those of observed satel- 
lite galaxies, which are thought to be associated with these 
subhaloes, has led to two apparent inconsistencies: (i) The 
radial number density distribution of subhaloes reveals a 
constant density core, whereas galaxies in clusters seem to 
follow an NFW profile, and (ii) whereas dark matter haloes 
of different masses look homologous when it comes to the 
properties of their dark matter subhaloes, a galaxy clus- 
ter looks very different from a galaxy sized system when it 
comes to their satellite galaxies (see Section 1 for references) . 
In this paper we have made use of the conditional 
luminosity function (CLF) to address these issues using 
data from the 2dFGRS. We constructed large Monte Carlo 
Markov Chains that allowed us to explore the full poste- 
rior distribution of the CLF parameter space. Using detailed 
mock galaxy redshift surveys (MGRSs) based on the CLF 
we have analyzed the radial distribution and abundances of 
host and satellite galaxies in the 2dFGRS. Our main con- 
clusions can be summarized as follows 

(i) The 2dFGRS is missing a significant fraction of galax- 
ies in close (projected) pairs (about 2 to 4 percent of all 
galaxies) due to the overlap and merging of galaxy images 
in the APM catalogue (see also Cole et al. 2001). 

(ii) Due to this close-pair incompleteness we can not 
put strong constraints on the radial distribution of satel- 
lite galaxies. After modelling the close-pair incompleteness 
in our MGRSs the data are consistent with a model in which 
the radial number density distribution of satellite galaxies 
follows that of the dark matter particles (i.e., a NFW dis- 
tribution), but we can not rule out alternatives with, for 
example, a constant number density core. 

(iii) Within the ACDM concordance cosmology with erg = 
0.9 we can not simulteneously match the abundances of host 
and satellite galaxies in the 2dFGRS. Matching the number 
of satellites per host requires exceptionally high mass-to- 
light ratios on cluster scales. This is in excellent agreement 
with previous findings based on the CLF formalism, but 
severely overpredicts the abundance of both host and satel- 
lite galaxies. 

(iv) Simultaneously matching the luminosity function, 
the clustering properties as function of luminosity, and the 
abundances of host and satellite galaxies seems to require a 
reduction of the power spectrum normalization to erg ~ 0.7. 
As shown by Yang et al. (2004a,b), this is also required in or- 
der to match the pairwise peculiar velocities in the 2dFGRS, 
the clustering power on small scales, and the multiplicity 
function of galaxy groups. 

(v) The CLF models for the ACDM concordance cosmol- 
ogy with erg = 0.7 that match the abundances of host and 
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satellite galaxies are extremely well constrained. They in- 
dicate that the average mass-to-light ratio of dark matter 
haloes reveals a pronounced minimum of 97 ± 3h (M/L)0 
at a halo mass of (7.2 ± 0.7) x W^°h~^ Mq. On cluster 
scales (M >, 10^*/i~^ M©) the average mass-to-light ratio is 
57O±4Ofe(M/L)0. 
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